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ABSTRACT 

We consider the role of advection in the two-temperature accreting corona with an 
underlying optically thick disc. The properties of coronal solutions depend significantly 
on the description of advection. Local parameterization of advection by a constant 
coeficient S replacing the radial derivatives lead to complex topology of solutions, 
similar to some extent to other advection-dominated accretion flow solutions. One, 
radiatively cooled branch exists for low accertion rates. For higher accretion rates two 
solutions exist in a broad range of radii: one is radiatively cooled and the other one 
is advection-dominated. With further increase of accretion rate the radial extensions 
of the two solutions shrink and no solutions are found above certain critical value. 
However, these trends change if the local parameterization of advection is replaced by 
proper radial derivatives computed iteratively from the model. Only one, radiatively 
cooled solution remains, and it exists even for high accretion rates. The advection- 
dominated branch disappears during the iteration process which means that a self- 
consistently described advection-dominated flow cannot exist in the presence of an 
underlying cold disc. 

Key words: galaxies: active - accretion, accretion discs - black hole physics - binaries 
- X-rays; stars - galaxies:Seyfert, quasars - X-rays. 



1 INTRODUCTION 

There are direct observational evidences that the gas sur- 
rounding central black holes in active galactic nuclei (AGN) 
and in galactic black holes (GBH) consists of two phases: 
colder, optically thick phase and hotter optically thin phase 
(for a review, see Mushotzky, Done & Pounds 1993; Tanaka 
& Lewin 1995; Madejski 1999). The emission of radiation is 
clearly powered by accretion onto the black hole. However, 
it is still not clear which of the two phases is responsible for 
accretion. 

In the classical disc/corona models (e.g. Liang & 
Price 1977; Bisnovatyi-Kogan & Blinnikov 1977; Haardt & 
Maraschi 1991) the accretion proceeds through the disc and 
the coronal gas does not contribute significantly to the angu- 
lar momentum transfer, presumably because of its magnetic 
coupling to the disc. 

In clumpy accretion flow models the accretion pro- 
ceeds predominantly through the cold clumps of gas (Collin- 
Souffrin et al. 1996, Czerny & Dumont 1998; Krolik 1998; 
see also Torricelli-Ciamponi & Courvoisier 1998). 

However, there are also models in which the hot gas car- 
ries most of the mass. The division of the flow into the two 



phases is predominantly radial in models based on advection 
dominated accretion flow (ADAF) solutions (e.g. Ichimaru 
1977; Abramowicz et al. 1995; Narayan, Kato & Honma 
1997; Esin, McClintock & Narayan 1997), in the Compton 
cooled solutions for ion tori (Shapiro, Lightman & Eardley 
1976; hereafter SLE), or in models considering both cool- 
ing mechanism (Bjornsson et al. 1996) . In those models the 
accretion proceeds through the cold accretion disc in the 
outer region and through the hot plasma region in the inner 
region. 

In this paper we discuss the model of an accreting 
corona. Initial formulation of the basic assumptions was 
given by Zycki, Collin-Souffrin & Czerny (1995) and the 
final formulation of the model was outlined by Witt, Czerny 
& Zycki (1997; hereafter WCZ). In this model the accretion 
proceeds both through a disc and a corona, in proportions 
determined by the model and varying with the distance from 
the black hole. The model was applied to Nova Muscae 1991 
(Czerny, Witt & Zycki 1996) and AGN (Czerny, Witt & 
Zycki 1997). 

In this paper we address the problem of advection in 
the corona. An optically thin accreting corona should show 
similar general behavior as a general optically thin flow not 
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accompanied by the disc, i.e. we might expect both advec- 
tion dominated solutions and radiatively cooled solutions. 
Advection was included in the model by WCZ and was found 
to be always negligible. Here we explain the nature of this 
phenomenon. 



2 MODEL 

2.1 Corona structure 

We assume that the corona itself accretes, i.e. that the en- 
ergy in the corona is due to the direct release of the grav- 
itational energy in the hot phase, without the necessity of 
having a mediator (e.g. magnetic field) between the disc and 
the corona. We further assume that the energy release can be 
described by the a viscosity prescription of Shakura & Sun- 
yaev (1973), with the energy generation rate proportional to 
the gas pressure (we neglect radiation pressure in the corona 
because the corona is optically thin). 

We assume a two-temperature plasma in the corona as 
in the classical paper of SLE, i.e. the ion temperature is dif- 
ferent from the electron temperature. The energy balance 
is computed assuming that the release of gravitational en- 
ergy heats the ions, the Coulomb coupling transfers this en- 
ergy to electrons and finally electrons cool down by inverse 
Compton process, with the disc emission acting as a source 
of soft radiation flux. The hot corona is radiatively coupled 
to the disc, as described by Haardt & Maraschi (1991). We 
assume isotropic emission within the corona (rj = 1/2) and 
for the (energy integrated) disc albedo we adopt a value 
a = 0.2. Compared to our previous paper (WCZ) we now 
use an accurate prescription for the Compton amplification 
factor, A (defined by F c = AF so ft), namely we compute 
A from Monte Carlo simulations of Comptonization. Our 
Monte Carlo code employs the method described by Pozd- 
nyakov, Sobol & Sunyaev (1983) and Gorecki & Wilczewski 
(1984). Assuming slab geometry (Thomson optical depth r es 
and electron temperature T c ) and the soft photons spectrum 
as a black body of temperature T s , we compute A on a, grid 
of T c , t cs , andT s and interpolate for values of interest at 
each radius. 

We neglect the dynamical term in the corona since it 
was shown to be relatively unimportant (WCZ). Therefore 
we can assume that the corona is in hydrostatic equilibrium 
and instead of solving original complex differential equations 
we adopt a set of simplified equations given in Appendix D 
of WCZ, with modifications concerning the advection and 
the Compton amplification factor. We also neglect the effect 
of the vertical outflow of the gas from the corona discussed 
by WCZ, i.e. we assume that the total accretion rate (the 
sum of accretion rates through the disc and the corona) is 
constant, independent of radius. 



2.2 Disc— corona transition 

The main feature of the model is that the division of the 
accretion flow into the hot and cold phases is not arbitrary. 
The thermal instability in an irradiated medium (Krolik, 
McKee & Tarter 1981) results in its spontaneous stratifica- 
tion into two stable phases. The criterium for the transition 



from cold to hot phase is formulated in terms of a specific 
value of the ionization parameter H which we define as 

(1) 

LJ gas 

where r/F c is the fraction of the coronal flux directed towards 
the disc and P gas is the coronal gas pressure (see also Kro- 
lik 1998). Following Begelman, McKee & Shields (1983) we 
adopt the following scaling, 

S = 0.65(T C /10 8 K)~ 3/2 (2) 

where T e is the electron temperature of the corona at a given 
radius. 

This description of the disc/corona transition can be 
easily argued for in a qualitative way. The essence of the 
transition is a change from Compton cooling mechanism 
to atomic cooling, including bremsstrahlung. Therefore, the 
transition zone corresponds to certain assumed contribution 
of bremsstrahlung to the total cooling. The criterion (Eq. ^|) 
applies accurately to the systems with either Compton heat- 
ing of the corona or a heating proportional to the density 
and fixes the bremsstrahlung contribution to 2/3 of the total 
cooling. Similar criterion formulated by Krolik (1998) gives 
the value of 3/7. Since the relative bremsstrahlung contri- 
bution decreases rapidly with the departure from the transi- 
tion zone into the corona, while the pressure remains roughly 
constant, we use this criterion as the basic criterion for pres- 
sure and we neglect bremsstrahlung as a cooling mechanism 
(see also Krolik 1998). 

This additional equation enables us to actually deter- 
mine the fraction of energy, /, which is liberated in the 
corona (N.B. in our previous paper WCZ we used £ = 1 — /). 
Since in our model / describes also the fraction of the mass 
accreting through the corona, we have 

M c (r) = f(r)M; M d (r) = [1 - f(r)]M (3) 

where Ma, M c and M are the disc, coronal and the total 
mass accretion rates. The total accretion rate does not de- 
pend on r if mass loss (e.g. through a wind as in WCZ) is 
neglected. 

The formulation of the model allows for computing the 
corona structure independently form the disc internal struc- 
ture. The disc/corona coupling is through surface pressure, 
Pg as , the coronal radiation flux, F c and the disc soft flux, 
F so f t uniquely specified by / and albedo. There is no need 
for subsequent iterations between the computations of the 
corona structure and the disc vertical structure as long as 
the local disc emission is well approximated by a blackbody 
(more generally: a thermal emission with a constant ratio 
of the colour and effective temperatures) and the albedo is 
fixed. When the local corona parameters are determined, 
they provide the surface boundary conditions for the equa- 
tions of the cold disc structure. The disc vertical structure 
can be solved if the viscous transfer within the disc is spec- 
ified (Rozariska et al. 1999). The imposed boundary condi- 
tions are automatically satisfied by the model, and the disc 
structure is determined uniquely. 

The radial variation of the relative proportion of the 
disc and coronal accretion flows may mean that either a frac- 
tion of the coronal mass cools and settles down on the disc 
surface thus joining the disc flow instead of falling radially 
into black hole, or a fraction of the disc flow evaporates from 
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the disc surface and joins the coronal flow. These changes 
are forced by the requirement of the hydrostatic and thermal 
equilibrium between the disc and the corona at each radius. 
However, the dynamics of this phenomenon is beyond the 
scope of our present model. 

2.3 Advection 

In the present paper we include the advection term in the 
corona in all computations. The advection is described as in 
WCZ, i.e. the energy balance takes the form 



F c = Fee 
and 



^rF adv = f(r)Mc 2 s S; 



dlnP 5dlnT; 



(4) 



(5) 



d In r 2 dlnr ' 
where F c is the energy flux generated in the corona, Fee is 
the Compton cooling of the corona by the soft disc photons 
and c a is the sound velocity in the corona. Equation |B| can 
be converted to give 



Ti 

T\r\r 



where T v j r is the virial temperature, 
_ _ GM m H 



k 



(6) 



(7) 



For the numerical solution of WCZ, 8 = 0.75 (see Appendix 
C in that paper). In general, however, 8 is a function of 
radius and should be computed consistently. In the next 
Section we will show solutions for a number of fixed val- 
ues of 8, while in Section ^ we will discuss solutions with S 
determined consistently through iterations. 

The parameters of the model are: the viscosity param- 
eter a in the corona and the dimensionless accretion rate rh 
measured in the Eddington units 

M 

M Edd = 3.52 __ [M Q /yr] (8) 

where M is the mass of the central black hole and we as- 
sumed the pseudo-Newtonian efficiency of accretion equal 
1/16 and pure hydrogen opacities. 

We show the results for the value of the viscosity pa- 
rameter a — 0.3 and the mass of the black hole M = WMq 
but we discuss the trends of solutions with these parameters 
varied. 

2.4 Spectra 

At each radius, the equations of the structure of the corona 
determine the electron temperature, the optical depth of 
the corona and the soft photon flux from the disc. The ef- 
fect of the Comptonization of the disc flux by the corona 
is calculated at each radius separately, using semi-analytical 
formulae of Czerny & Zbyszewska (1991). We neglect here 
the anisotropy of the Compton scattering which should be 
taken into account in detailed spectral modeling (e.g. Pouta- 
nen & Svensson 1996, Haardt, Maraschi & Ghisellini 1997). 
However, very accurate computation of the spectra is not 
the main goal of the present paper. 

The final disc spectrum is computed by integration over 
the disc surface assuming an inclination angle equal zero 



(i.e. top view). This integration procedure is an essential el- 
ement of our model since the corona properties are radius- 
dependent: outer parts of the corona are predominantly re- 
sponsible for the high energy extension of the spectrum and 
its hard X-ray slope while the inner, cooler parts of the 
corona mostly influence the soft X-ray range by producing 
moderately comptonized disc component. 

All computations are done for a non-rotating black hole 
and the relativistic corrections are neglected. 



3 CORONA PROPERTIES FOR CONSTANT S 

In this Section we show results obtained assuming a value 
of the advection parameter 8 (Eq. ||) constant with radius. 



3.1 The relative strength of the corona 

The most important prediction of the model is a strong ra- 
dial dependence of the strength of the corona. This depen- 
dence changes qualitatively with the accretion rate, mainly 
due to the presence of advection. Examples of radial depen- 
dences of the fraction of energy generated in the corona, 
/(r), are shown in Figure [l] for a number of values of 8 and 
rh. 

Spatial extent of the corona is always finite, indepen- 
dently of rh and 8. The corona covers only an inner part of 
the disc from a certain outer radius i? max inwards. 

At i? max all the energy is liberated in the corona for low 
accretion rates, if 8 (and for any rh when 8 = 0, i.e. if 
there is no advection). Consequently, at r = i? max the ac- 
cretion flow proceeds entirely through the corona. At larger 
radii no corona solutions of our equations exist since the 
Compton cooling provided by the disc is too large, under 
the adopted assumptions about the corona structure. There 
is therefore a strong and discontinuous change of accretion 
flow structure at i? max . For r > fl max all the accretion pro- 
ceeds through the disc whilst at r = -R max the accretion 
proceeds through the corona, with the cold disc heated only 
by X-ray illumination. We envision that the change occurs 
through rapid heating and evaporation of the disc surface 
which takes place in order to maintain the thermal balance 
condition, but the detailed dynamics of this process is be- 
yond the scope of our present model. Closer in, the relative 
strength of the corona decreases and, consequently, the rel- 
ative fraction of the disc accretion increases. 

The dependence of the radial extension of the corona, 
Rmnx, for 8 = 0.75 on the accretion rate is shown in Figure 
^| The size of the corona is considerable, covering the disc 
up to ~ 130J?schw for accretion rate approaching 0.04 MEdd, 
but 7? ma x decreases significantly for smaller accretion rates, 
down to about 10J?schw for sources radiating at 0.1 per 
cent of the Eddington luminosity. For accretion rates below 
4 x 10 -4 A^Edd the corona ceases to exist, as the Compton 
cooling is too str ong while heating becomes inefficient, (see 
also Section 5.3). 



Although the relative fraction of energy dissipated in 
the corona increases with radius (up to 7i max ) the actual 
energy flux decreases outwards. We can see that from the 
simple analytical solutions given in Appendix C of WCZ. 
Since / increases with the radius r approximately as / oc 



r 3/8 



the X-ray flux Fx(r) decreases with radius: Fx( 
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Figure 1. The fraction of the energy dissipated in the corona, /, as a function of radius for viscosity parameter a = 0.3, mass of the 
black hole IOMq, and four values of the advection parameter 8 (Equation (5 = 0, 0.2, 0.5 and 0.75. Labels inside panels refer to 
accretion rate. 



r ' . This may be important for detailed computations 
of relativistic smearing of the X-ray reprocessed component. 
However, we do not discuss this spectral component in the 
present paper. 

A second branch of solutions appears when the ac- 
cretion rate rh approaches a certain critical value, m a< jv 
On this branch the cooling is dominated by advection and 
rh a dv is a function of 5: m a d v = 0.192, 0.069, 0.044 for 



8 — 0.2, 0.5, 0.75, respectively (Figure jjj). The exact topol- 
ogy of the two solutions is a rather complex function of <5 and 
rh. For rh just above m a d v the advective solution coexists 
with the radiative one in a range of radii. The two branches 
cross for somewhat higher rh and then separate, creating two 
spatial regions where both solutions can exist, with an inter- 
mediate range of r, where no solution is possible. The outer 
ring then shrinks rapidly as rh increases and disappears, the 
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Figure 2. The dependence of the extension of the corona on 
the accretion rate for 8 = 0.75. Light shades mark the Compton 
cooling dominated solutions and the dark shades mark the region 
of both advection dominated and cooling dominated solutions. 
The viscosity parameter a = 0.3 and the black hole mass is IOMq. 
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Figure 3. Coronal solutions shown in coronal accretion rate, m c , 
versus surface density in the corona, E, computed at 5 fischw 
for a = 0.3. Open circles show advection-dominated solution 
and solid squares show the Compton cooling dominated solution 
for 8 = 0.75. Note that the advective branch terminates when 
F a dv/Fc = 0.9, i.e. strongly advective solutions do not exist in our 
model. Solid triangles show the solution with <5(r) computed con- 
sistently through iterations (Section ^) . Only one solution branch 
exists in this case, and it extends into fairly high rh c (correspond- 



ing to total M > M Edd ). Black hole mass M = 10 M 



©■ 



inner one does the same but more slowly. Mathematically, 
the solutions continue to (unphysical) / > 1 region creating 
closed loops. 

For comparison we also plot in Figure |l| the radial de- 
pendence of the relative strength of corona obtained when 
the advection term was neglected (5 = 0). Obviously, in this 
case only one branch of solutions appears. It is similar to 
the Compton-cooled branches of solutions with 8 > for 
low accretion rates, rh £ 0.01, i.e. the advection is then 
only a small correction to the energy balance. For larger 
rh the 5 = solution is quite different from the advection- 
corrected solutions: the corona extends to much larger radii 
and no qualitative change with increasing rh is seen. It is 
clearly the advection term which is responsible for the com- 
plex topology of solutions seen in Figure [j]. 

The merging of the radiatively cooled solution with the 
advectively cooled solution has been previously found by 
Chen et al. (1995) and Bjornsson et al. (1996) for an op- 
tically thin flow, not accompanied by any disc. Therefore, 
we also constructed an analogous log rh — log E diagram for 
our coronal solution at 5i?schw (Figure |^). We plot coro- 
nal surface density E c versus the coronal accretion rate, 
rh c = rhf(r) where the value of factor /(r) depends on 
the solution branch. 

The two solution branches for a given total rh pro- 
duce two points on this plot, since / is different on the 
two branches. The two solutions merge for rh = 0.09 (at 
our assumed radius 5i?schw, ct = 0.3 and 8 = 0.75), result- 
ing in a continuous curve with both the fraction of energy 
carried by advection and the fraction of energy dissipated 
within the corona increasing with frh. The uppermost point 
is characterized by / = 1, i.e. the whole energy is generated 
within the corona, but advection transports ~ 90 per cent 
of the energy rather than 100 per cent. In our model the 
fully advective branch does not develop, due to the pres- 
ence of the disc providing soft photons for Compton cool- 
ing. Mathematically, the fully-advective branch would ap- 



pear for F BO f t = 0, which would require unphysical condition 
/ = 1/[1 — 77(1 — a)] > 1. Thus the necessary presence of the 
soft photons from the disc suppresses in our model the fully 
advective branch, presen t in other optically thin solutions, 
e.g. ADAFs (see Section 5.1). 

The existence of the advection dominated branch of so- 
lutions was not found by WCZ, and the reason for that is 
given in Section 3.2.2. 

The viscosity parameter a has a strong influence on the 
existence and properties of the corona. For a lower than 
our assumed 0.3 the solutions are more limited in rh and 
radius, the more so the higher the advection coefficient 8. For 
example, for a = 0.03 and 8 = 0.2 only very spatially limited 
solutions exist. For the same a but 5 = 0.5 no solutions exit. 



3.2 Corona properties 

3.2.1 Advection 

In Figure ^ we show the ratio of the energy flux advected 
inwards with the coronal flow to the total dissipated flux 
for 8 = 0.75. For low accretion rates only the radiatively 
cooled solution exists and advection is not important. How- 
ever, when rh approaches rh = 0.1, up to 40 per cent of the 
flux on this branch is carried by advection. This fraction 
depends on the radius. 

For high accretion rates (rh ^ 0.04) the second solution 
appears. This solution is cooled mostly by advection and it 
is similar to ADAF solutions. Since in that case most of the 
energy (more than 90 per cent) is liberated in the corona 
the accretion proceeds mostly through the corona itself. 

3.2.2 Ion temperature and the geometry 0} the corona 

The ion temperature decreases almost inversely with radius 
(sec Figure 5 in WCZ and formulae in Appendix C in that 
paper). The pressure scale height of the corona, defined as 
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Figure 4. The radial dependencies of the ratio of the advected energy flux to the total flux generated in the corona (upper left panel), 
optical depth (upper right panel), electron temperature (lower left panel), and ion temperature (lower right panel), for viscosity parameter 
a = 0.3, mass of the black hole 10A/q, and S = 0.75. The thin straight line in the lower right panel shows the virial temperature. The 
dashed line shows advection-dominated solution whilst the solid line shows the solution cooled predominantly by Compton cooling. 
Labels inside panels refer to accretion rate. 



\GMmn ) 



(9) 



increases almost linearly with radius and the ratio Hp / R is 
almost constant. The coronal accretion flow resembles ac- 
tually a spherical accretion, similarly to the case of pure 



ADAF flows. In such flows the ion temperature is of the or- 
der of the virial temperature and geometrical thickness of 
the flow is of the order of r. Sound velocity is close to the 
Keplerian velocity and is proportional to the radial velocity 
, where the proportionality coefficient is given by viscosity 
parameter a. However, our coronal solutions are generally 
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cooler and Hp/R ratio, although constant, is not equal 1. 
Nevertheless, the dependence of the radial velocity on radius 
is quite similar to the case of ADAF. In Figure ^| we plot 
the ratio of the radial to sound velocity, computed from the 
formula: 



Vr = 



fM 
47ri?E c 



(10) 



We see, that far from the marginally stable orbit v r /v s ra- 
tio is of order of a and is almost constant throughout the 
disc. The flow is then moderately subsonic, depending on the 
value of the viscosity parameter. Close to the marginally sta- 
ble orbit our coronal flow, as well as ADAF solution, become 
transonic and continue as a free fall onto a black hole. 

In this paper we use a simplified description of the verti- 
cal hydrostatic equilibrium and we have to check afterwards 
whether the solution can actually be in hydrostatic equilib- 
rium. The approximate criterion is that the ion temperature 
should be smaller than the local virial temperature. We see 
from Figure ^ that 71 is usually lower than T v i r on the 
radiative branch, but 71 > T v i r on the advective one. The 
same problem refers to the cor ona height as a function of ra- 
dius. Since Hp / R = ^T\/T v i r , the ratio Hp / R can be larger 
than 1, i.e. the corona can be (very) geometrically thick. The 
reason for 71 exceeding Tvir can be seen from Equation ^ 
Ti = (Tvir/5) x (-Fkdv/^c), i.e. 71 can approach and exceed 
T v i r if advection is dominant. 

The super-virial ion temperature is the reason why 
advection-dominated solutions were not found by WCZ. In 
that paper the vertical structure was calculated much more 
carefully, assuming the hydrostatic equilibrium at the basis 
of the corona and allowing for the transonic vertical outflow 
from the corona. These solutions automatically prohibited 
the violation of the hydrostatic equilibrium at the basis of 
the corona. In that case the set of solutions for increasing 
accretion rates simply ended up as soon as the ion temper- 
ature reached the virial temperature (see Figure 2 of WCZ 
and Section 3.1 and 4.1 therein). 

The same problem was noted by Narayan & Yi (1994) 
in the case of pure ADAF solutions and expressed as the 
problem of the Bernoulli constant being positive for ADAF. 
Therefore, the model with a hot medium being in hydro- 
static equilibrium does not offer correct solution beyond cer- 
tain value of accretion rate. 



3.2.3 Electron temperature and the optical depth 

Since the density in the corona decreases outwards the effi- 
ciency of Coulomb interaction between the ions and elec- 
trons decreases as well. However, both the disc and the 
corona bolometric luminosities go down rapidly with the ra- 
dius. Therefore the electron temperature, T e , rises outwards 
and the ratio of Ti/T e decreases outwards. The highest value 
of T e is of the order of 1.5 x 10 9 K, and it depends only weakly 
on m and the viscosity parameter, a. Such a universal value 
is an interesting property of our model for lower accretion 
rates. At higher accretion rates, however, the outer radius 
of the corona contracts rapidly due to the advection and the 
maximum corona temperature also rapidly drops down with 
an increase of the accretion rate. 

The optical depth of the corona is practically indepen- 
dent of radius, and very weakly dependent on other param- 
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Figure 5. The ratio of the radial velocity to the sound velocity in 
the corona as a function of radius. The approximate constancy of 
v T /v B for r 10i?schw means that the geometry of the accretion 
resembles spherical accretion. The viscosity parameter a = 0.3, 
mass of the black hole IOA/q and two values of the accretion rate 
m: m = 0.01 (continuous line) and m = 0.001 (dashed line). 
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Figure 6. Coefficient S(r) computed from the solution for m = 
0.2, and initial So = 0.2 (equation |i ) . Solid squares represent the 
radiative branch, while the circles show the advective branch, af- 
ter the first iteration (Si). On the advective branch the computed 
S is systematically larger than the initial one which leads to the 
disappearance of this branch, when solution with consistent S(r) 
is calculated. Solid triangles (Sao) show S(r) for the converged 
solution, where only one branch remained. 



eters; it is always between 0.1 and 0.2 (see Appendix C in 
WCZ). 



4 SOLUTIONS WITH CONSISTENT S(r) 
4.1 Corona structure 

The assumption that the advection parameter S is constant 
as a function of radius is not correct in general. As can 
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be seen from equation |5| 8 is determined by radial deriva- 
tives of the ion temperature and pressure, hence it can be a 
function of radius. Since the topology of solutions depends 
rather sensitively on S (Figure we can expect significant 
changes of the topology if we consistently compute the 5(r) 
dependence. 

As a first step, we show in Figure ^| the coefficient 5 
computed numerically from the solution obtained for ini- 
tial constant So = 0.2 and rh — 0.2 (labelled Si) For these 
parameters two solution branches exist: with radiative and 
advective cooling dominant (cf. Figure |l]) . On both branches 
the computed 8 has a strong radial dependence. However, 
as Figure |l] shows again, the radiative branch is present for 
any 5, therefore we do not expect significant changes to its 
character, when solution with proper S(r) is computed. The 
same is not true for the advective branch: the higher the S, 
the narrower the range of this branch's existence in the r-rh 
plane. Since the computed 5 is now significantly larger than 
the initial one, we can expect a decrease of the importance 
of the advective branch. 

Proper self-consistent solutions describing the flow can 
in principle be obtained either by solving radial differential 
equations containing explicitly the advection term expressed 
as derivatives given by equation |B|, or by an iterative proce- 
dure correcting the distribution of S(r) at each iteration step. 
First method was successfully applied by Chen, Abramow- 
icz & Lasota (1997) and Narayan, Kato & Honma (1997) to 
obtain global solutions for ADAF flow. We used the second 
method, iterating solving the corona structure equations for 
a given <5(r) and computing the next approximation to S(r). 
Our algorithm is very similar to the method used by Chen 
(1995): in order to find solution at a given radius r, we solve 
the equations at r and two auxiliary radii, r — Ar and r + Ar 
(assuming <5(r± Ar) = 5(r), but the solution is insensitive to 
this particular condition) . We then compute the corrected S 
from equation |^ and compute the corrected structure. When 
convergence is achieved, we proceed to the next radius. 

Iterating the solution for the advection-dominated 
branch quickly leads to its disappearance. The radiative 
branch shrinks somewhat (i? max decreases), especially for 
rh such that the advected fraction close to -R max was sub- 
stantial in the non-iterated solution. In other words, only 
the solutions with rather low advective cooling survive. 

Figure ^ (points labelled 5oo) shows an example of the 
iterated 5(r) dependence. The iterated S(r) is positive for 
larger radii, where advection is a cooling process, but it 
changes sign for smaller radii, as advection becomes a locally 
heating process. The same trend was obtained for optically 
thick discs calculated taking into account advection, depar- 
ture from the Keplerian rotation and the transonic character 
of the flow close to the marginally stable orbit (Muchotrzeb 
& Paczyhski 1982; Abramowicz et al. 1988). The heating role 
of advection increases dramatically close to the marginally 
stable orbit since the energy dissipation there approaches 
zero. 

In Figure |^ we show examples of the converged solu- 
tions of the corona structure as functions of rh and a. The 
corona is strongest at is outer edge, although / is not al- 
ways 1 at r = -Rmax- Similarly to solutions with constant 
<5, the electron temperature is ~ 100 keV while the optical 
thickness is ~ 0.15. Advection is never important as a cool- 
ing process. At small radii, r ^ 20 -Rschw, 5 changes sign, 



so the advective flux contributes to heating. In this region 
the solution exist for rh ^ 4 x 10 -4 , and up to rh c ~ 1, 
i.e. the Eddington luminosity in the corona. The topology 
of solutions represented on the rh-E diagram change as well 
(Figure ^): the solution forms a single branch only. 

The solutions are sensitive to the value of the viscosity 
parameter a. For low a the solutions are generally rather 
limited spatially. Moreover, for a ^ 0.1 the ion temperature 
strongly exceeds the virial temperature, so the solutions can 
hardly be considered physically acceptable. 

The value of the mass of the central black hole influences 
the results only very slightly. Solutions for M = 10 8 Mq 
practically overlap those for M — 10 Mq. This is to be ex- 
pected (see e.g. Bjornsson & Svensson 1992 and references 
therein), since the only dependence on the central mass in 
our model is through the temperature of the soft flux, which 
can affect the Compton amplification factor, but the depen- 
dence A(To) is very weak for steep (soft) spectra. 

4.2 Radiation spectra 

Spectra predicted by the model are rather soft and domi- 
nated by the disc emission. Examples are plotted in Figure 
^| for parameters characteristic for AGN and GBH. Gener- 
ally, the power law components are harder for lower rh. 

For accretion rates of order of 0.05 AfEdd, typically ex- 
pected in Seyfert galaxies, the amount of hard X-ray emis- 
sion predicted by the model is negligible and the high en- 
ergy index far too steep, so the present model does not offer 
a promising description. On the other hand, the original 
model of WCZ (with neglected advection) reproduced well 
the typical properties of AGN (see Czerny et al. 1997 for 
radio quiet quasars and Seyfert galaxies, and Kuraszkiewicz 
et al. 1999 for NLS1 galaxies; see also the corona parameters 
determined for MCG+8-11-11 by Grandi et al. 1998). 

For the Galactic black holes the general trends are simi- 
lar to those for AGN. For very low accretion rates (rh Ss 0.01) 
the model predicts considerable hard X-ray emission from 
the corona, although the spectra are somewhat steeper than 
those resulting from the original model by WCZ (Janiuk & 
Czerny 1999). Very steep spectra are predicted for m > 0.1, 
i.e. corresponding to the high or very high state of stellar 
black hole systems, while V ~ 2 in the observed spectra of 
e.g. soft X-ray transients (Zycki, Done & Smith 1998). 



5 DISCUSSION 

In this paper we considered the accreting corona model. 
Main features of the model and differences from other hot 
accretion flow solutions are, firstly, the presence of a cool, 
optically thick disc supplying soft photons for Comptoniza- 
tion in the corona, and secondly, the condition for (vertical) 
stratification of the flow into the two phases due to thermal 
instability. Since in this model the accreting coronal plasma 
is hot, optically thin and geometrically thick, it is necessary 
to include the radial energy transport (advection), similarly 
to the optically thin ADAF solutions and optically thick, 
slim discs. 

The structure of solutions crucially depends on whether 
the advective flux is solved for consistently - i.e. the coeffi- 
cient 5 (Eq. ^) is a function of radius - or S(r) is assumed 
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Figure 7. Examples of the corona structure computed with consistent solution for the advection coefficient 5(r). Two values of the 
viscosity coefficient a are shown, with two values of rh for each a, as labelled. The solutions were computed for 10 Mq, but the 
dependence on the BH mass is negligible. 



constant, although there are certain properties of the solu- 
tions independent of the ^-prescription. 

For a constant S, at a given radius (smaller than a cer- 
tain maximum value) there is either one solution - radia- 
tively cooled through Compton process - or two solutions, 
one radiatively cooled and the other advection-dominated. 
For radii larger than the maximum one no coronal solutions 
exist. Whenever two solutions are possible, there is always 
a critical value of the accretion rate for which the two solu- 
tions merge and no solutions are found for higher accretion 
rates. 

Similar effect is seen when studying the radial struc- 
ture of coronal solutions. At very low accretion rates (below 
rh ~ 4 x 10 -4 ) there are no coronal solutions. At larger rh 
the radiatively cooled solution appears, which covers an in- 
ner part of the disc. This solution was previously found by 
WCZ. The fraction of the disc covered by the corona in- 
creases with the accretion rate but the fraction of radiation 
emitted by the corona decreases. At even higher rh (above 
rh w 0.04, depending on 5) a second, advection-dominated 
solution emerges, if the corona structure equations allow for 
the ion temperature to be higher than the virial temper- 
ature. The two solutions merge at the outer edge of the 



corona. With further increase of the accretion rate, the re- 
gion covered by the corona contracts rapidly, with no corona 
present for rh <; 0.1 (for the viscosity parameter a = 0.3). 
Smaller a leads to solution merger for even lower accretion 
rates, as found previously by e.g. Abramowicz et al. (1995) 
and Bjornsson et al. (1996). 

The existence of the advection-dominated solutions was 
automatically excluded in the original formulation of the 
model by WCZ - instead, they simply observed the dis- 
appearance of the single, radiation cooled solution with an 
increase of the accretion rate due to the ion temperature 
reaching the virial value. 

When the <5(r) function is solved for, we obtain only 
one solution branch. It is cooled by radiation i.e. with ad- 
vective cooling negligible. In fact, advection changes sign for 
r ^ 20 -Rschw he. it acts as a heating process. The proper- 
ties of this branch are similar to those of the radiatively 
cooled branch obtained for a constant 8: the solution ex- 
ist for rh > 4 x 10 -4 , and the corona is strongest at its 
outer edge. For radii such that advection is a heating pro- 
cess there, this solution can exist up to rather high rh, even 
formally exceeding MEdd in the corona. It disappears only 
when the temperature of the disc flux reaches the (decreas- 
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Figure 8. The radiation spectra for M = 10 Mq (left panel) and M = 10 s Mq (right panel). Two values of the accretion rate are 
plotted, as labelled, and the viscosity parameter a = 0.3. Corresponding solutions of the corona structure are plotted in Figure H. 



ing) electron temperature in the corona, but it is strongly 
super- virial already for lower m. 



5.1 Comparison with other hot, optically thin 
solutions 

Our solutions show general similarity to other hot, optically 
thin disc solutions (SLE; Abramowicz et al. 1995; Chen et al. 
1995; Bjornsson et al. 1996; Narayan & Yi 1995; Zdziarski 
1998; see Kato et al. 1998 for general discussion). There 
are also, however, certain important differences, as a direct 
consequence of the assumptions specific to our model. 

When plotted on the rh-S diagram (Figure ^), our 
solutions for constant 5 form two branches, merging for cer- 
tain maximum m. Both the existence of the two branches 
and their merging was found previously. However, as op- 
posed to 'conventional' ADAF solutions, strongly advection- 
dominated branch does not appear in our model, due to 
substantial Compton cooling. For the same reason our solu- 
tions disappear for lower viscosity, a ^ 0.03, while no such 
effect is observed in the above mentioned solutions, where 
the adopted cooling mechanism can usually be made suffi- 
ciently inefficient (see below). 

We obtain an increase of the maximum allowed coronal 
rh, when we solve for the advective flux (i.e. 5(r)). This has 
also been found previously by Chen (1995). Again, however, 
the advection-dominated branch does not appear in our so- 
lution, and the solution disappears for low viscosity. None 
of these effects has appeared in the Chen (1995) work. 

Hot discs cooled by comptonization of external soft pho- 
tons were considered by Zdziarski (1998), as a generalization 
of the SLE solution. However, the source of soft photons 
was not specified in that work i.e. the soft flux could im- 
plicitly be assumed to be small, if required. Therefore the 
general properties of the solutions found were in close agree- 



ment with the previous solutions where bremsstrahlung was 
usually assumed as the radiative cooling process. For ex- 
ample, the strongly advective branch is present in Zdziarski 
(1998) solution, in spite of the efficient radiative cooling. The 
Compton parameter y = AkT e /m D (? r is ~ 1 in that work 
(it is ~ 0.1 in our solutions), which indeed requires very low 
Fsoft in order not to suppress the advectively-cooled branch. 

Where an ADAF-like disc is assumed to co-exist with 
a cold, optically thick disc (e.g. Esin et al. 1997), the tran- 
sition radius is an adjustable parameter. In our work the 
additional equation (Eq. |l|) resulting from the thermal insta- 
bility condition closes the structure equations, thus enabling 
us to compute the transition radius. 



5.2 Super-virial ion temperature 

The ion temperature in an optically thin ADAF flow tends to 
be larger than the local virial temperature and the Bernoulli 
constant for such a flow is positive. Therefore such solutions 
are not possible without some kind of outflow (e.g. Narayan 
et al. 1997, Blandford & Begelman 1999). 

The same problem affects the advection dominated 
branch of our coronal solutions and the iterated solutions 
obtained for low a. Such a corona violates the assumption 
of the hydrostatic equilibrium. Coronal gas with ion tem- 
perature exceeding the local virial temperature cannot flow 
in, as assumed in our model or ADAF solutions. If strong 
outflow indeed developed, no accretion would take place, 
switching off the energy source. A moderate transonic out- 
flow from the corona surface perpendicular to the disc does 
not provide a solution. It was already included in the original 
formulation of the model by WCZ and it did not prevent the 
disappearance of coronal solutions for accretion rates higher 
than m ~ 0.1. A magnetic wind can provide a solution, 
provided that the outflow launched at certain radius carries 
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away more angular momentum than what is locally neces- 
sary for accretion to proceed (Blandford & Begelman 1999) . 
However, at this stage such solutions are rather arbitrary. 



5.3 Corona formation at its outer edge 

Our model does not provide a mechanism for disc evapo- 
ration. Instead, it allows to check whether the corona, if 
formed, can exist in hydrostatic and thermal equilibrium. 
The outer edge of our corona is therefore the maximum ra- 
dius at which those conditions are satisfied. Since we did not 
consider two-dimensional flow, our transition from a bare 
disc to a coronal solution is sharp. The transition would be- 
come smooth if the radial conduction were included (e.g. 
Honma 1996), but this would require two-dimensional com- 
putations, since the radial thickness of the transition zone is 
expected to be same as the vertical thickness (see discussion 
by Dullemond 1999). The coexistence of the bare disc and 
the disc/corona at i? ma x does not contradict the analysis of 
Dullemond & Turolla (1999) since the coronal part is not 
strongly advection-dominated, with approximately half of 
the energy (or less) carried by advection and the remaining 
energy radiated away locally, as in SLE. 

The mechanism leading to corona formation is still un- 
specified. It may be related either to disc instabilities, or 
magnetic phenomena. The coronal solutions are generally 
within the standard disc radiation pressure instability zone 
(e.g. Janiuk & Czerny 1999), but this does not seem a strong 
argument in favour of the first possibility. 

The transition radius, i? ma x, in our models is gener- 
ally rather smaller than the outer radius of the ADAF flow 
in typical applications, which is assumed ~ 10 4 -Rschw (e.g. 
Esin et al. 1997). Since our description of the flow applies 
in principle to the outer part of an ADAF solution, where 
the hot flow and the cold disc flow are assumed to over- 
lap (Esin et al. 1997), it may mean that large ADAF radii 
would also be difficult to achieve if the disc/corona coupling 
was correctly included. However, our result should be rather 
treated as an indication of a possible problem and not a def- 
inite answer. The model depends rather sensitively on the 
adopted description of the disc/corona coupling, expressed 
through Equation ^, as was shown for a non-advective 
corona by Janiuk & Czerny 2000. The simple change from 
the bremsstrahlung contribution from 2/3 (our assumption) 
to 3/7 (Krolik 1998) would only change the results quantita- 
tively by radially expanding the region of coronal solutions 
but lowering the optical depth of the corona. However, better 
description of the transition region with disc evaporation in- 
cluded may change the results more significantly. It may also 
provide an explanation of the formation of ADAF part. Un- 
fortunately, although some preliminary, partial results are 
available (e.g. Meyer & Meyer-Hofmeister 1994, Dullemond 
1999, Rozariska & Czerny 1999), they are still not directly 
applicable to the models and further research is needed. It 
may also be true that the radial extension of the overlapping 
region will appear to be small, as suggested by Dullemond 
(1999), thus further complicating the prediction of location 
of the disc/ADAF transition. 



6 CONCLUSIONS 

• Advection-dominated branch of coronal solutions does 
not represent a physically acceptable description of the flow. 

• Accreting corona solutions are pre- 
dominantly Compton-cooled and, for small radii, exist for 
all accretion rates larger than ~4x 10 _4 MBdd- 

• Spectral slopes predicted by accreting corona models 
for AGN and GBH are too steep in comparison with obser- 
vations so the disruption of the innermost part of the disc 
or magnetically driven outflow seem to be required. 
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